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lL , 

" I suggest that the great body of knowledge gained over the past 10 years about sim- 

<y) i pie spin- 1/2 quantum quantum antiferromagnets points to a connection between 

cuprate superconductivity and the strong interactions. The underlying physical 
idea, which I admit to be highly speculative, is that the phase diagram of such 
a magnet consists of competing ordered phases regulated by a nearby quantum 
critical point. Exactly at this critical point the low-lying elementary excitations of 
i-^j | the magnet are gauge fields and particles with fractional quantum numbers anal- 

(-H , ogous to the spinon and holon excitations found in spin chains. An arbitrarily 

small distance away, however, these bind at low energy scales to make the familiar 
collective modes of the ordered states into which one renormalizes. Vestiges of 
these "parts" of the collective modes may be seen in conventional materials and 
models in high-energy spectroscopy and inconsistencies in sum rules exactly the 
way quarks are seen in particle physics. 
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^ ■ Figure 1: The physical behavior of the strong interactions generated spontaneously in simple 

antiferromagnets. 

k>( \ The premise of this article is illustrated in Fig. 1. I wish to argue that 

; i ■ there are physically identifiable objects in simple Heisenberg antiferromagnets 

which behave like U(l) quarks and could conceivably be apt analogues of them. 
These are not the elementary excitations of the system in most cases but rather 
objects out of which the elementary excitations are built. It is my current belief 
that the quark-like objects and the gauge fields through which they interact are 
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the true elementary excitations at some nearby quantum critical point, but I 
shall mostly sidestep this issue and concentrate on the physical meaningfulncss 
of the particles in the commonly-studied cases. 

The antiferromagnets in question are described by the t-J Hamiltonian 

H = P G H P G , (1) 

where 

p o = I[{ 1 - c h c n c n^} (2) 

j 

is the Gutzwiller projector and 

n o= E{-*£4 c ^+i s ;- s 4 > ( 3 ) 

the sum < j,k > being over near-neighbor pairs of a lattice, with each pair 
counted twice to maintain hcrmiticity. When the dimension of the lattice is 2 or 
greater, the phase diagram of this model is complex and beyond our means to 
compute reliably. When the dimension of the lattice is 1, however, there is an 
exact solution at the supersymmetric point J = 2i, the ground state of which 
is a nondegenerate singlet, i.e. has no order, and the elementary excitations of 
which arc spin-1/2, charge-0 particles known as "spinons" and spin-0, charge-1 
particles known as "holons" . These are the quarkJike objects of the problem. 
Their dispersion relations are plotted in Fig. 2 a. The existence of spinons 
and holons is intimately connected with the lack of order in the ground state, 
and is thus common in 1 dimension, where continuous symmetry breaking is 
impossible, but uncommon in higher dimension, where order seems to occur 
almost always. But if the higher-dimensional ground state is forced to be 
disordered by means of a variational ansatz, which is equivalent to adding a 
small long-range interaction to the Hamiltonian to destabilize the order, then 
spinons and holons makes sense and we obtain in 2-d the dispersion relations^ 

E „ = l. 6 J v / cos 2( fe)+cos 2 ( ^ ) (4) 

E hoion = ± 2t^cos 2 (q x ) + cos^q y ) (5) 



2 



-TT o n r A s r 




Figure 2: Left: Spinon and holon dispersion relations obtained by Bares, Blatter, and Ogata 
from the Bethe ansatz solution of the supersymmetric spin chain. Right: 2-d spinon and 
holon dispersion relations given by Eqs. (5) and (6). Inset: 2-d brillouin zone. 

plotted in Fig. 2. These are also implicit in the U(l) gauge thearj: descriptions 
of the t-J model based on the commensurate flux saddle point Era. I wish now 
to establish that these particles have physical meaning at intermediate energy 
scales even when the system is allowed to order, giving rise to forces that 
bind them at low energy scales into the well-known excitations of the ordered 
phases. 

In Fig. 3 I show the optical conductivity 

<r*x(u) = ~ V | < a\j x \0 > \ 2 S(hiu -E a + E ) (6) 

a 
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Figure 3: Left: Optical conductivity computed by exact diagonalization by Moreo and 
Dagotto on a 4 X 4 cluster for the case of S = 1/16 and J/t = 0.4. The arrow at lu = 
indicates the Drude oscillator. Right: Total oscillator strength — <T> /Nt defined by Eq. 
(9) and its Drude contribution computed by exact diagonalization on a 4 X 4 cluster. The 
symbols □,+, and O correspond to J/t = 0.1,0.4, and 1.0, respectively. The solid curve is 
a plot of Eq. (11). The dashed curve is a guide to the eye. 



computed by Dagotto and MoreoB for a single hole in a 4 x 4 cluster, which is 
representative of such calculations. Here \a> indicates an exact eigenstate of 
energy E a and j x is the electric current operator 

3x = Pg t| YsY, | c \° Ck ° ~ c k* c i° f Pg ' ( ? ) 

j (X > 

where k denotes the near neighbor of j in the x-direction. This calculation 
provides evidence for the existence of the holon and measures the size of its 
mass. The /-sum rule 

( w ^ <omo> 

a xx {u)dh) = — (8) 

is plotted versus doping in the lower part of the figure, as is its "Drude" 
contribution. The width and shape of the later cannot be computed accurately 
but its integrated area can. Both sum rules are straight lines at low doping, 
the slope of which does not depend on J. This is the behavior of a doped 
semiconductor. From the usual expression 
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we find a mass of 



Jo 



n n 

a xx (u})duj = - — n , (9) 
2 m 



h 2 

m = Q.71— , (10) 
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where b is the bond length. This compares favorably with the 1 / v2 in these 
same units obtained from the curvature of Eq. (5) near its minimum and the 
0.54 obtained in 1-d from the Bethe solution. The fact that the "Drude" weight 
is always about half the total indicates that this particle is the carrier. The 
full sum rule 

<T>= -2.6 NtS (11) 

also agrees with Eq. (5) in equaling the — v8£ per hole associated with the 
holon band minimum. This number has no connection to the mass in general, 
and is thus an additional constraint on the band structure. 

Further evidence for the existence of the holon may be found in the electron 
propagator in the limit of small J/t. Following the notation of Eq. (7), we 
define the electron propagator by 



GV(£)=£{ 



< a| C t g |Q > | 2 | < a\c qa \0 > | 2 



E - E a + E Q + ij] E + E a - E - %r\ 



(12) 



where 

N 



7= J2exp{iq- r 3 )c, ja . (13) 



j 

In Fig. 4 I show the imaginary part of this function at half-filling calculated by 
the exact diagonalization method for J/t — 0.0 and 0.2 by DagottoQ. In either 
case the spectrum is a broad continuum about 6t wide with a pronounced 
dip in the center and a weight that moves from low to high energy as the 
momentum is advanced from T to M. In the J — > limit the broad continuum 
may be ascribed to the decay of the injected hole into spinon-holon pair in the 
limit that the spinon is very heavy, for then the spectrum should be the holon 
density of states 

^) = T^E I"'" S(E - E h q t° n ) dq x dq y (14) 

weighted by a decay matrix element. The density of states computed from Eq. 
(5) is plotted in Fig. 4. Any reasonable model will give the motion of the 
weight since 

< Olc^lO >=i (15) 
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Figure 5: Left: Quasiparticle dispersion relation calculated by Liu and Manousakis using 
spin-wave perturbation theory for the case of J/t = 0.2. Right: Quasiparticle bandwidth 
W/t calculated by Poilblanc using exact diagonalization on clusters of various sizes. The 
dashed line is a plot of Eq. (17). 



<0\c\ a H c 9CT |0>=<0|7i|0> -It 

at half-filling, where 1 and 2 denote near-neighbor sites. 

The J/t = 0.2 curves also have a peak at low binding energy which is 
the quasiparticle of the magnetic insulator. In Fig. 5 I show the dispersion 
relation of this quasiparticle found numerically by a number of authors u. It 
has a deep minimum at E and an overall bandwidth W, the difference between 
the maximum and minimum of the dispersion relation, that does not depend 
nn t. This width, measured in multiples of t, is plotted against J/t in Fig. 5 
Q. From the slope of the line one obtains 



<Si -S 2 > 



cos(q x ) + cos(q y ) V (16) 



W = 2.2J , (17) 

or 1.6^/2 J, which is the spinon bandwidth given by Eq. (4). The prefac- 
tor 1.6 in Eq. (4) has the physical significance of a magnetic stiffness. It 
causes the spinon velocity at £ to be the spin-wave velocity of the ordered 
antiferromagnetl 

«. = 1.6^ . (18) 

The quasiparticle peak is accompanied by scattering resonances. These 
cannot be seen in Fig. 4 because the sample is too small, but they may be 
seen clearly in Fig. 6, which is the spectral pfunction at £ for J/t = 0.2, 
calculated using spin wave perturbation theorp. The quasiparticle peak and 



7 



J/t 




-10 12 3 4 

E/t 



0.1 0.2 0.3 0.4 0.5 
J/t 



Figure 6: Left: Spectral density at S calculated by spin wave perturbation theory by Liu 
and Manousakis for the case of J/t = 0.1 in the limit of large sample size. Right: Energies of 
quasiparticlc (I) and first two string resonances (II and III) as a function of J/t. The dashed 
lines are plots of Eq. (19). 



the first two resonances are labeled by roman numerals. Their energies are 
plotted as a function of J/t in Fig. 6. The lines through the data points 
represent the formula 



E n /t = -3.28+ (J/t) 2/3 x 



2.03; n = 1 
5.46; n = 2 
7.81; n = 3 



(19) 



These energies are exactly the spectrum expected a light particle in orbit about 
a heavy one, provided the attractive force between the two is a string, i.e. 
V(r) ~ \r\. More precisely, the Hamiltonian 



W = -^V 2 + 2.2 J\'-\ 
2m b 



-3.28i 



(20) 



where m is the mass derived by the conductivity sum rule and given explicitly 
by Eq. (10), has energy eigenvalues given by Eq. (19) except for substitution 
(2.63, 5.54, 7.81) -» (2.03, 5.46, 7.81). 

These facts have the following physical interpretation. The quasiparticle 
is a bound state of a spinon and the holon analogous to the hydrogen atom. 
Its band structure tracks that of the spinon because the spinon is "heavier" 
than the holon in the sense of having a narrower band. The optical sum rule, 
by contrast, is sensitive to the light particle, and thus measures the holon 
properties. The same thing is true in hydrogen, where the acceleration mass 
is dominated by the proton but the optical properties are dominated by the 
electron. The potential binding these particles together is a string at low 
doping, which means that they can never separate and do not exist as separate 
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entities in this limit, but already at a doping of one hole in a 4 x 4 lattice, 
or S = 1/16, something occurs to allow the string to break and the holon to 
ionize off to become a free carrier. 

The t-J Hamiltonian is formally equivalent to the Lagrangiantla 



TV 



ihdt + ' 



fia + b 



ihdt + 1 



£ \-j\Xjk\ 2 + oXjJ 

<j,k> 



J2fj<yfka ~%b]b k 



J 



(21) 



where fj a and bj are fictitious fermion and boson operators on site j in terms 
of which the electron is written 



'"0° ~ J 3 aU j ' (22) 
t>j is a Lagrange multiplier which when integrated out forces the constraint 



£/£/*. + >>■>> 



i 



(23) 



and Xjk is a Hubbard-Stratonovich variable. This is a U(l) gauge theory to 
the extent that \jk niay be approximated as having a fixed length, for then 
the phase functions as a vector potential along the bond < j,k >, the scalar 
potential on site j being <f>j . This turns out to be a bad approximation for this 
particular Lagrangian, but we can imagine adiabatically transforming it into 
one for which \xjk\ is fixed and for which a small Maxwell term 



C Max = -<J ^2 XjkXkiX£mXmj + -j ^2 ^ d tXjk + {(p 3 -<f>k)X]k\ 2 > (24) 

<j,k,im> <j,k> ' 

is added as a regulator. Then the classical saddle point has magnetic flux 7r 
per plaquette, the f- and b-particles become free particles with the relations of 
Eqs. (4) and (5), although with different coefficients, and we obtain conven- 
tional lattice QED with doubled fermions. The limit relevant to the low-doping 
numerical work is g — > oo, which is strongly confining. Thus the string forces 
may be associated with confinement in strongly-coupled QED, the antiferro- 
magnetic order in limit may be associated with the chiral symmetry breaking 
known to accompany confinement in this problem, and the spin wave, which 
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is both the Goldstone of the broken symmetry and a bound pair of spinons, 
may be associated with the pion. 

The correct appearance of a string force in the antiferromagnetically or- 
dered phase suggests that the unbinding of the quasiparticle seen in Fig. 3 
might indicate a first-order transition to a superconducting phase correspond- 
ing to the coulombic phase of the gauge theory. The magnetic order is known 
to disappear at about 5 = 0.05, which is consistent with deconfinement by 
5 = 1/16. However it is only a suggestion, for the above Lagrangian is less 
accurate and more difficult to solve than the spin Hamiltonian from which it 
was derived, and all the major ordering questions for the former are still un- 
resolved. It should be viewed not as a computational tool but as means for 
understanding how the physics of the strong interactions might materialize in 
a quantum antiferromagnet without being postulated. 

I wish to express special thanks to E. Dagotto for providing me his unpub- 
lished J — > spectral functions and to A. M. Tikofsky for numerous helpful 
discussions. This work was supported primarily by the NSF under grant No. 
DMR-9421888. Additional support was provided by the Center for Materials 
Research at Stanford University and by NASA Collaborative Agreement NCC 
2-794. 
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